Chiral particle separation by a non-chiral micro-lattice 
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We conceived a model experiment for a continuous separation strategy of chiral molecules (enan- 
tiomers) without the need of any chiral selector structure or derivatization agents: Micro-particles 
that only differ by their chirality are shown to migrate along different directions when driven by 
a steady fluid flow through a square lattice of cylindrical posts. In accordance with our numerical 
predictions, the transport directions of the enantiomers depend very sensitively on the orientation 
of the lattice relatively to the fluid flow. 
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Many building blocks of life (e.g. a ligand) are or- 
ganic molecules of one specific chirality which interact 
very differently with other molecule species (e.g. a recep- 
tor) due to their handedness (enantiomers): While one 
species may act as a drug, its chiral partner may be inert 
or even toxic [1]. Methods for separating enantiomers are 
thus of great importance in fields ranging from molecular 
biology and organic chemistry to pharmacy and agricul- 
ture. Up to date, all those methods require suitable chiral 
derivatization agents or some other kind of chiral selec- 
tor material or structure [H . Their main disadvantage is 
that they have to be developed and precisely predesigned 
anew for every given chiral analyte. 

Here, we put forward an alternative approach which 
does not require any such selector with a built in struc- 
tural chirality, but instead exploits a dynamical chiral 
symmetry breaking. The proposed model experiment 
from Fig. [T] builds on previous separation strategies for 
non-chiral particles, as conceived theoretically e.g. in 
Refs. [3|,|| and realized experimentally e.g. by Austin, 
Sturm, and coworkers [s^. Conceptually, our approach 
may be traced back to de Gennes, pointing out that 
small (but still macroscopic) objects of opposite chiral- 
ity generically must exhibit different transport directions 
when sliding (or rolling) down a solid plane or floating 
on a liquid surface under the action of a constant force 
^ . More specifically, chiral particles which are propelled 
in a shear flow have been explored in several theoret- 
ical 0] and experimental (sl-llH studies, demonstrating 
that the right- and left-handed particles experience a lift 
force of opposite sign in the vorticity direction of the 
flow. For particles of one specific chirality, the result- 
ing deflection from the main flow direction and its de- 
pendence on the shear rate was reported by Markino et 
al. for millimeter-sized chiral objects [10], and by Mar- 
cos et al. for micrometer-sized helical bacteria [llj. A 
genuine separation of right- and left-handed, centimeter- 
sized crystals was demonstrated by Howard [sj . Another 
main theoretical concept established by Kostur et al. ex- 
ploits the interplay of thermal noise and of a flow field 
with spatially variable vorticity to separate stylized chi- 



ral "molecules" |l2|. In a next step, Eichhorn theoreti- 
cally predicted that right- and left-handed particles are 
pumped with different velocities through a microfluidic 
channel with asymmetric walls [13i] . 




FIG. 1: Schematic top- view (x-^-plane) of the microfluidic ex- 
periment with fluid- and particle-reservoirs at its ends. The 
magnified inset shows an optical micrograph image of the cen- 
tral, microstructured area (square lattice of cylindrical post) 
and two particles of opposite chirality. A steady fluid flow in 
x-direction is generated by a pneumatic pump (not shown), 
maintaining a constant pressure difference between the reser- 
voirs. The angle a denotes the inclination of the square lattice 
relatively to the fluid flow. 

With our present work we further pursue yet another 
theoretical proposal [l^, namely to employ periodic po- 
tentials for separating chiral particles. Our experimental 
realization of this concept is schematically illustrated in 
Fig. [TJ featuring a central microstructured square lat- 
tice array extending over 1 cm in x-direction and 1 mm 
in ^/-direction, which consists of cylindrical posts with 
3.3 ± 0.2 lira radius, 20.5 ± 0.2 /im lattice constant, and 
6±0.2/im height (z-direction). The entire device is filled 
with distilled water, containing L- and F-shaped "enan- 
tiomers" in sufficiently high dilution so that their inter- 
action can be safely neglected. Their long and short axes 
are 14.7 ± 0.3 /im and 10.2 ± 0.3 /im respectively, and their 
height (z-direction) is 3.5±0.3/im. Hence, they are small 
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enough to edge their way through the lattice and large 
enough to prevent their swapping of chirality by rotation 
about their long axis. For a more detailed description 
of the particle and microstructure fabrication we refer to 
[iBj . The experiment is actuated by a pneumatic pump, 
maintaining a constant pressure difference of p ~ 15mbar 
between the reservoirs (see Fig. [1]), resulting in a steady 
fluid flow of the order of 100/im/s. During the experi- 
ment, a small part of the structured area is observed with 
an inverted microscope (Zeiss Axiovert 200) and the par- 
ticles are tracked by digital video microscopy. 

As emphasized in [l^, there does not exist any sym- 
metry argument of why the L-and F-shaped enantiomers 
should travel with identical mean velocities v through the 
square lattice when driven by a steady flow. Following 
de Gennes @, we thus can conclude that the two ve- 
locities will indeed be different unless we can name any 
a priori reasons to the contrary. One important such 
reason identified in [14] are locking phenomena in the 
absence of thermal fluctuations Since thermal noise 
effects turn out to be rather weak in our present experi- 
ment, the velocities of the two enantiomers are expected 
to be strictly speaking different but still practically indis- 
tinguishable under many conditions. To avoid such un- 
wanted cases and rather optimize the experimental set-up 
apriori, we performed extensive numerical simulations 
jig . In particular, we numerically explored many differ- 
ent (but still experimentally feasible) post diameters, lat- 
tice constants, particle dimensions, pressure differences, 
and inclination angels a. In the following, we immedi- 
ately focus on the resulting "best case prediction" , from 
which we derived the set-up specified in Fig. [T] and in 
the previous paragraph. 
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FIG. 2: Deflection angles i9r and of the F- and L-shaped 
particles versus lattice inclination a, predicted by means of 
numerical simulations for parameters corresponding to the 
experimental set-up from Fig. [T] The error bars arise by 
averaging over many initial conditions and realizations of the 
noise in combination with a long but finite temporal average. 



The quantity of foremost interest is the time- averaged 
particle velocity and in particular the deflection an- 
gle of this velocity from the x-axis. Those numerically 
predicted deflection angles for both the L- and F-shaped 
particles versus the inclination angle a of the lattice are 
presented in Fig. [2l For a < 9°, one of the above men- 
tioned locking effects occurs: Both particle species es- 
sentially follow the (1,0) lattice direction ('^^l = '^^r = <^), 
apart from extremely small deviations (not resolvable in 
Fig. [2]) due to thermal noise. In contrast, for a > 9° quite 
appreciable differences i^i, — i^r between the two transport 
directions arise. Beyond a = 16°, the error bars of 'dr 
and beyond a = 21° also those of of i^i, dramatically in- 
crease in our simulations. The reason is that, depending 
on the initial conditions, numerical solutions with differ- 
ent time-averaged velocities v may coexist in the absence 
of thermal noise. Since thermal noise is rather weak in 
our system, the velocities still exhibit a large variance 
within our quite long but still finite simulation times. To 
avoid such an undesirably slow convergence, we confine 
ourselves to a < 16°. For symmetry reasons it is further- 
more sufficient to consider a > 0. 

According to Fig. [2l the largest splitting angles 
— between the L- and F-shaped particles are ex- 
pected around a = 15°. A slightly weaker splitting, but 
remarkably enough with an opposite sign of '^^l — '^^r, is 
expected around a = 12.5°. Therefore, we will focus 
in the experiments on these two most promising values 
a = 15° and a = 12.5° predicted by the theory [l7|. 

As already mentioned, we performed similar simula- 
tions as in Fig. [2] for many other lattice and particle 
parameters. In order to maximise I'^^l — '^r|, the param- 
eters from Fig. [2] seem to be (practically) optimal, but 
significant — '^^rl could still be obtained for a wide 
range of substantially different parameters. 
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FIG. 3: Experimentally observed trajectories of L- and F- 
shaped particles for the same set-up as in Fig. [1] with a = 
15°. The total length of each trajectory depends on where 
the particle happened to enter and exit the visual field of the 
optical microcope. 

Fig. [3] illustrates how the L- and F-shaped particles 
actually move in the experiment for a = 15°. The main 
point of Fig. [3] is also the main point of our present 
paper: In accordance with our numerical predictions, the 
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two chiral particle species can indeed be separated from 
each other in our experimental microfluidic device with 
extremely high reliability. In fact, not a single particle 
with a "wrong" sign of i9 was observed. 
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FIG. 4: Left: Stroboscopic images by superimposing several 
video microscopy snapshots for one F- or one L-shaped par- 
ticle from Fig. [3] (a = 15°). The motion proceeds from left 
to right at time steps dt. Right: Corresponding numerical 
simulations. 



such transitions between rows of posts (see also the blue 
trajectories in Fig. [3]), but at a substantial smaller rate 
than in the simulations. In addition, for some of the L- 
shaped particles we also observed a different type of mo- 
tion in the experiment, illustrated by Fig. [5^, without a 
comparable theoretical counterpart. (While the particles 
in Figs. HJd and [5^ move into identical directions, their 
speeds (note the different time-steps dt) and orientations 
differ.) Quantitatively, we thus obtained an experimental 
deflection angle of = 12.1° ± 3.4°. 
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FIG. 5: (a): Same as in Fig. Od (a = 15°), but for a different 
video sequence of an L-shaped particle and a different time 
step dt =0.72 s. (b): Same as in Fig. [6^ (a = 12.5°), 
but for a different video sequence of a F-shaped particle (the 
snapshots are blurred due to the rapid particle motion). 

Fig. m provides a more detailed picture of how the par- 
ticles travel through the square lattice. The stroboscopic 
snapshots at time steps dt =0.48 s in Fig. demonstrate 
a very regular, time-periodic motion of the F-shaped par- 
ticles in the experiment. This behavior is extremely well 
reproduced by the simulations in Fig. Hh. The reduced 
time step of dt =0.24 s in Fig. \^ in addition reveals a 
quite pronounced rotational see-saw motion (the corre- 
sponding experimental snapshots are even more blurred 
than in Fig. [Sb and are therefore not shown). Due to 
thermal noise and imperfect spatial periodicity, also very 
rare deviations from the behavior in Fig. 3^ are observed 
(see also red trajectories in Fig. [3]), resulting in a mean 
experimental deflection angle of 'dr = —10.7° ± 0.8°. 

Turning to the L-shaped particles. Fig. |4]3 evidences 
an extremely "straight" motion along one and the same 
inclined row of posts. The simulation in Fig. |4]i also 
shows short episodes of such a straight motion, but ad- 
ditionally every now and then the particle bumps with 
its "short leg" against a post, turns around, collides with 
another post, and finally continues along the next lower 
row of posts with another short episode of straight mo- 
tion. In fact, the complete experimental video, of which 
only a small part is covered by Fig. IDd, also exhibits 



FIG. 6: Left: Same as in Fig. H but for a = 12.5°. 

Similarly, Fig. [6] illustrates the experimentally and 
numerically observed particle motion for the second of 
the above predicted "most promising" deflection angles, 
namely a = 12.5°. In contrast to a = 15°, now the 
agreement between experiment and theory is excellent 
for the L-shaped particles, but somewhat worse for the 
F-shaped particles. Moreover, for the F-shaped particles 
there is again a second type of motion in the experiment 
without a theoretical counterpart, see Fig. [5]3. Alterna- 
tively, the motion of the F-shaped particles predicted by 
the theory may be viewed as a compromise between the 
two experimentally observed types of motion. 

We emphasize once again that the main purpose of 
our numerical simulations was to forecast a promising 
layout for our subsequent experimental realization of a 
microfluidic chiral selector device. In view of Fig. [3l 
our theoretical approach was clearly very successful in 
this respect. Accordingly, a further improvement of the 
simulated model [l6| is beyond the scope of our present 
work. Yet, since the quantitative agreement between sim- 
ulations and experiments in Figs. [2]l6]is clearly not fully 
perfect, it is worthwhile to point out some plausible rea- 
sons. First of all, it is obvious, both from the simula- 
tions and the experiments, that quite small changes of 
a may lead to quite notable changes of the particle ve- 
locity with respect to direction, modulus, and the co- 
existence (or not) of more than one "type of motion". 
Numerically, a comparably sensitive dependence of the 
particle velocities was also observed for most of the other 
model parameters. One therefore expects that even quite 
small imperfections of the theoretical model may easily 
result in very significant errors of the predicted particle 
motion. Indeed, our model [l6| involves several approx- 
imations of this kind, most notably: (i) Along the z- 
direction, the microfluidic device from Fig. [T]is bounded 
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by a "bottom-" and a "top-plane", both of which are, 
however, not perfectly planar but rather exhibit notable 
imperfections/roughness/defects. Every now and then, a 
traveling particle "bumps" into one of them (due to ther- 
mal diffusion in 2:-direction etc.) with the result of a quite 
notable, sudden "deflection" (even visible under the mi- 
croscope), (ii) The experimental post array deviates from 
a perfect square lattice and the single posts are not at all 
of perfect cyhndrical shape, (iii) While the force exerted 
by the fluid flow on the particles is estimated in the the- 
oretical model by means of the fluid flow in the absence 
of the particles, in reality there will be quite notable cor- 
rections due to hydrodynamic interactions between par- 
ticles and microstructure. In principle, it is not difficult 
to account for all those effects (i)-(iii) in a theoretical 
model dynamics, but its numerical simulation would re- 
quire (among others) to solve the full, three-dimensional 
hydrodynamics at every time step, which goes far beyond 
of what is feasible in practice. On the other hand, the 
qualitative and quantitative agreement between theory 
and experiment in Figs. [2J6] is still remarkably good in 
view of all those approximations (i)-(iii). 

Conclusion: The main point of our paper was an ex- 
perimental proof-of-principle that chiral particles can be 
separated without employing any kind of chiral selec- 
tor structure but rather by exploiting a dynamical chiral 
symmetry breaking. In order to predict a suitable design 
of the experimental set-up, a simple theoretical model 
dynamics has been introduced and numerically solved. 



Any further improvement of the various approximations 
of this model itself are beyond the scope of our present 
work. As illustrated by Fig. [3l not a single particle in the 
experiment exhibited a deflection angle with a "wrong" 
sign, i.e. the chiral separation fidelity was extremely 
high. A more systematic exploration and optimization 
of the proposed microfluidic device, e.g., with respect to 
some specific technological applications or a quantitative 
comparison with other methodologies will certainly be 
an important issue in the future, but goes beyond the 
objectives of our present study. 

Besides the chiral separation per se, one of our main 
findings was that the motion of the chiral particles 
through the periodic microstructure depends very sen- 
sitively on the specific values of many dynamical and 
structural system parameters. On the one hand, the par- 
ticle motion is thus very difficult to predict qualitatively 
by means of simple intuitive arguments and with high 
quantitative accuracy by means of numerical simulations. 
On the other hand, this sensitivity underlines the practi- 
cal efficiency and variability of the experimental concept. 
The most promising future direction is now to downscale 
the generic principles behind our present approach into 
the realm of chiral nano-particles and molecules. 
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